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The early history of the Vogt-Russell theorem is retraced following 
its route starting at the realization of a correlation between mass and 
luminosity of binary and pulsating stars, through the embossing 
of this observation into a theorem, and finally to the emerging first 
signs of its failure to serve as a theorem in the strict mathematical 
sense of the word. 


Introduction 

Astrophysics is not an exact science in the sense of mathematics. 
On the observational side, time and again the Universe is consid¬ 
erably more creative than human imagination; on the theoretical 
side, the equations that are derived to model the intricate processes 
and phenomena in astrophysics usually lack the simple symme¬ 
tries or properties to attract mathematicians' attention for in-depth 
analyses of their properties. Therefore, proclaiming a theorem in 
astrophysics - by astrophysicists - was, is, and likely will remain a 
daring undertaking. 


THE THEOREM IN STELLAR PHYSICS, the topic of this margina¬ 
lia, is known by the name of Vogt-Russell (VR) theorem: It claims, 
roughly speaking, that the structure of any star is uniquely de¬ 
termined by its mass and its chemical composition alone. Today 
we know that the theorem does not hold in its original strict for¬ 
mulation; nonetheless, the VR theorem continues to enjoy some 
popularity and pops up in almost all courses on stellar structure 
and evolution, in textbooks, and even recently in research papers. 

^ The following exposition retraces the history of the VR theorem 
during the (semi-) analytic era of stellar astrophysics. The first 
statement of what not much later became the theorem appeared 
essentially en passant in a research note of Heinrich Vogt where 
he generalized Eddington's mass-luminosity relation. The 'proof' 
by Henry Norris Russell, which advanced the original claim to a 
theorem, was enshrined later in Chandrasekhar's seminal reference 
work. An Introduction to the Study of Stellar Structure. After that 
the VR theorem lived an apparently quiet and unquestioned life 
to the end of the 1950s when first counterexamples, albeit rather 
un-astrophysical ones, were put forth. The literature of the early 
1960s contains more cases of by then more naturalistic star models 
that seemed to violate the VR theorem. The following exposition 


covers the history up to the publication the monograph of Cox & 


Giuli (1968^ where the authors devoted a short chapter to the VR 


* e.g. |Carroll & Ostlie|j20i^ or 
|Melis et al.|j2oi4| 
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theorem and gave an impression of the thinking on this matter at 
that time. 

A forthcoming second part of the history of the VR theorem will 
deal with a revived interest in the topic in the 1970s, an epoch when 
stellar evolution theory had turned into a computation-intensive 
enterprise with ever more complex input physics and ever more 
complicated resulting star structures. 


First formulations 


The 1910s and 1920s were the years when astrophysicists started to 
understand the stars as long-living self-gravitating fluid spheres. 
Even though the source of energy was nof yet identified, the 
thermo-mechanical structure of stars was already modeled math¬ 
ematically. In accordance with this intellectual achievement, the 
body of history-related literature collecting, scrutinizing the con¬ 
tributions, and cross-linking the players in the field is as huge as 
authoritative. For those interested in the topic, Hufbauer ( |20o6| , 
Gingerich ( 1995I, Cowhng] ([i966[|, and references therein serve as 


fertile starting points. 

For this marginalia it is sufficient to realize that in 1924 Arthur 
Eddington published a paper with the title On the relation between 
the masses and the luminosities of the stars ( |Eddington|i924) where 
he used his analytical thermo-mechanical star models to fit the 
observed correlation of luminosities and masses of stars (see Ap¬ 
pendix A for a discussion). Eddington collected the data of 46 stars 
belonging to various kinds of binary stars for which masses and 
absolute magnitudes were reported; all of them (plus five pulsating 
variable stars) obeyed a remarkably smooth relation. 

In December 1925, Heinrich Vogt, an astronomer at Heidelberg's 
Kdnigstuhl Observatory, submitted a short theoretical paper to 
the Astronomische Nachrichten ( |Vogt|i9^ in which he generalized 
Eddington's analytical star models by introducing spatially variable 
forms of the mass-absorphon coefficient and of the energy gener¬ 
ation rate. In the last paragraph of the research note, Vogt men¬ 
tioned laconically, what he later referred to as the Eindeutigkeitssatz 
(uniqueness theorem): 


[...] Wir miissen annehmen, dafidie mittlere Dichte, die effektive Temperatur 
und die absolute Leuchtkraft eines Sternes nur von seiner Gesamtmasse 
abhdngen [...] 


Vogt allowed, additionally, for a small variation of the magni¬ 
tudes of the global stellar quantities at fixed mass because they 
might be influenced by the nature of the stellar material, i.e. by the 
star's chemical composition. In other words, Vogt claimed that the 
global stellar quantities L* (M, }i), Tgff (M, }i), p{M, }i) are uniquely 
determined by the value of stellar mass and the star's composition 
alone. As it seems, the uniqueness claim was not particularly im¬ 
portant to Vogt: About one year after the first note on the subject, in 



Figure 1: Heinrich Vogt (1890 - 1968) 
after [Bohrmann, A.|{1968) 


Translation: We must assume that the 
mean density, the effective temperature 
and the absolute luminosity of a star 
depend on its total mass only. 


The chemical composition, which 
is assumed to be homogeneous in 
the stars rmder consideration, be 
quantified by its mean molecular 
weight p. 
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a long review of the theory of stellar structure and evolution, Vogt 
( |i927^ did not even touch the aspect of uniqueness of the solutions 
to the equations describing the structure of the stars. Another year 
later, when Vogt ( [i928[ l devoted a more extensive paper parhcu- 
larly to the mass - luminosity law to emphasize once again that a 
form very close to that of Eddington can be recovered even if the 
internal structure of the stars differs from what Eddington assumed 
And again, the paper of 1928 does not mention the Eindeutigkeits- 
property. Only in 1930, in the longest treatise ( |Vogt|i93^ concern¬ 
ing the relation of mass, luminosity, and effective temperature of 
the stars - essentially an attempt to understand the distribution 
of the stars in the Hertzsprung-Russell (HR) diagram - did Vogt 
refer briefly and informally to the dependence of the global stellar 
quantities on mass and chemical composition; however, without 
mentioning his 1926 paper. Even after the 1930 publication, one 
wonders how much importance Vogt actually attributed to the 
uniqueness conjecture and if he realized its consequences. 

On the other side of the Atlantic, Henry Norris Russell studied 
the physical basis of stellar evolution since the earliest days of 
modern astrophysics (cf. DeVorkrn 20oo| . In the mid 1920s, he 
and colleagues at Princeton Observatory overhauled the textbook 
Astronomy - A revision of Young's Manual of Astronomy, updating 
it also with the latest ideas and results from the thriving field of 
stellar astronomy ( Russell et al.|i92^ . In the second volume, in 
Section 975, the authors state: 


[...] It is found that a star of given mass and composition will usually he in 
equilibrium for only one value of the radius, and hence for definite values of 
the luminosity and surface temperature. For stars of different masses these 
values will be different, but so long as the composition is the same, all the 
stars of a given luminosity will have to be of some one definite size, surface 
temperature, and spectral type. [...] 

The above statement was made without further elucidation and 
without going into any technical details. In any case, the unique¬ 
ness claim of Russell et al. appears more deliberately formulated 
and in parhcular physically more coherent than Vogt's statement 
from a year earlier. In the same Section 975 of the textbook, the 
uniqueness property was then applied to the interpretation of the 
structure of what later became the HR diagram and of Eddington's 
mass - luminosity diagram. Any scatter to the observed distribu¬ 
tion of stars in these diagrams that goes beyond the observational 
uncertainties was interpreted by the authors to mean that [...] gener¬ 
ation of heat [...] is different in different stars. The stars cannot therefore 
all contain the same proportion of'active matter' [...], in other words [...] 
they [the stars] must differ in composition. Hence, in the textbook of 
Russell et al. ( 1927), the uniqueness statement of the structure of the 


stars is of auxiliary use only, namely as a supporting argument in 
the interpretation of the distribution of the stars in the HR diagram. 
The stars were thought to populate the HR plane as a function of 


At the time of Vogt's article, the 
color/spectral-type - magnitude 
diagram was typically referred to as 
the Russell diagram. The designation 
'Hertzsprung-RusselT diagram took 
over in the literature only in the 
second half of the 1930s. 



Figure 2: Henry Norris Russell (1877 - 
1957) after|DeVorkin| (2013) 
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mass and of varying chemical composition. The mass was believed 
to diminish with age so that stars evolve across the HR plane from 
high to low mass during their life. At the end of Section 975 of 
Astronomy, where contributors to the content of the stellar-evolution 
discussion, including the uniqueness conjecture, were referred to, 
only Russell himself and Eddington appeared - Vogt's research 
note was not mentioned. 


The path to theorem 


In a review, which discussed the state of the theory of the consti¬ 
tution of the stars, Russell ( [i93ia| clarified that Vogt and, indepen¬ 
dently, he himself formulated a few years earlier a theorem on the 
uniqueness of stellar structure. When referring to Vogt, Russell 
cited the extensive paper of Vogt ( |i93o| on the nature of the correla¬ 
tions of global stellar observables and their relation to their internal 
structure rather than the short note of 1926. One volume of the 
MNRAS later, RusseU ( |i93ib) eventually set the records straight: 

In an addendum, he reported that Vogt brought to his attention the 
proper reference confarning the first sfafement of the uniqueness 
claim, which goes back to the year 1926. In any case, as early as 
1931, Russell regarded the uniqueness statement as a theorem (of 
deep insighf and as a confribufion fo fhe field of sfellar asfrophysics 
which he considered to have remained much undervalued in the 
community). Not much later, RusseU| ( [i933| | reiterated his opinion 
in a non-technical survey on stellar astrophysics in general, and 
stellar structure and evolution in particular; he referred again fo fhe 
Vogt theorem ^ as the [...] most important general proposition regarding 
stellar constitution [...]. As of fhen, Russell apparently attributed to 
the uniqueness theorem much more importance than Vogt ever did 
in any of his wrifing; this might be connected with Russell's interest 
to thoroughly interpret the distribution of sfars across the spectral 
class - luminosity diagram, on which he actively worked with the 
goal to pack all the stars into a coherent story in the framework of 
the theory of stellar evolution as it stood then (cf. Gingerich|i995 
Hufbauer|20o6). 


Russell, in contrast to Vogt, resorted already early on to a math¬ 
ematical argumentahon and proofed the uniqueness theorem (cf. 
Russell|i93ia) : The stellar structure problem was considered as 
the solution to a system of four differenhal equations with dis¬ 
tributed boundary conditions. Since the number of boundary 
condifions was counfed fo be fhree, compared wifh four differenfial 
equations, Russell concluded after some meandering that unique 
one-parameter sequences of solufions must exist for a prescribed 
chemical composifion, and that without loss of generality, the star's 
mass can be chosen as this parameter. 

Under the spell of fhe success of quanfum fheory fhe young Dan¬ 
ish astrophysicist Bengt Strdmgren 3 set out, in the early 1930s, to 
improve the understanding of the HR diagram adopting hydrogen- 


’ in part 2, p. 417 


In Section 10, Russell sets the stage 
by proclaiming self-confidently: 

[...] The proof is simple. [...] 


3 cf. |Rebsdorf|(2oo3) 
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rich stellar models with the use of more elaborate microphysics. 

In long review paper, Stromgren ( 1937) presented a comprehen¬ 
sive view of his understanding of the theory of the stellar interiors 
and of sfellar evolution. A whole section 4 of the exposition was 
dedicated to the uniqueness conjecture of Vogt and Russell. Strdm- 
gren referred to it as the [...] Satz 5 von Vogt und Russell [...} (p. 477) 
and he also outlined its proof, following the line of argumentation 
of Russell ( i93ia^ . Not much later and apparently influenced by 
Strdmgrens work, the authoritative monograph Stellar Structure of 
[Chandrasekhar 1 1939^ consolidated the theorem status of what, at 
best, should still have been considered the conjecture of Vogt and 
Russell. In Section Vll.i of his book, Chandrasekhar contemplated, 
nonetheless, circumstances which could void the VR theorem. He 
h5rpothesized material fimctions as sources of trouble; e.g. a nuclear 
energy generahon rate, which does not depend on the local values 
only of p and T. Macroscopic counterexamples, however, were yet 
beyond the intellectual horizon. 

(19531 ^ studied homology transformahons of the stellar 


Kurth 


structure equations and their properties (in a remarkably modern 
formulation). In this context, Kurth formulated an aggravated ver¬ 
sion of the VR theorem in that he progressed from if to iff: In the 
framework of homologous, chemically homogeneous star models 
in complete equilibrium 7 he concluded that the stars have nearly 
the same internal structure iff they have nearly the same mass and 
nearly the same chemical composition. Only a few sentences after 
this conclusion, Kurth cautioned the reader that [...] Beiviesen ist 
nichts, es handelt sich nur um Plausibilitdtsbetrachtungen [...], all state¬ 
ments were indeed clearly declared as plausibility considerahons; 
in particular, he explicitly assumed that solutions to the stellar 
structure equations exist. Finally, Kurth also distinguished between 
the pure stellar structure problem (i.e. stationary solutions to the 
structure equations), for which the VR theorem was formulated, 
and the full stellar-evolufion problem which, by its very nature, 
is a time-dependent problem. Therefore, he warned that any ex¬ 
trapolation of conclusions from the application of the original VR 
theorem to the realistic stellar-evolution problem has to be treated 
with utmost caution. 

Only a few years later, Odgers 11957' ^ reported his attempts 
to construct homologous series of chemically homogeneous star 
models in complete equilibrium. The aim was to derive a simple an¬ 
alytical formulahon of the mass - luminosity law of main-sequence 
stars. Already early on in the paper, Odgers criticized the mathe¬ 
matical assumptions that entered the uniqueness statement of the 
VR theorem and he even deduced an explicit homologous series 
which violated it. Because the inferred energy-generation law in 
the conflicting model series was unphysical in the stellar context 
the counterexample to the VR theorem was considered at best of 
formal interest. Even though Odgers' paper was an actual mathe¬ 
matical blow for the VR theorem, the paper made no impact in the 


■t Section II. 16 


5 Satz in the mathematical sense of 
theorem. 


^ Kurth, Rudolf, *1917 in Germany, 
mathematician; doctorate (1948) and 
habilitation (1951) at the University of 
Berne, Switzerland, on topics of stellar 
dynamics. Later he pursued an aca¬ 
demic career in England and the USA; 
he was a prolific writer of books on 
epistemology, philosophy, mathematics 
in physics and astronomy. 

^ Complete equilibrium in the sense of 
thermal and hydrostatic equilibrium. 


Translation: Nothing is proved, these 
are all plausibility considerations only 


* Odgers, Graham ]., *1921 in Australia; 
he spent his life as an astronomer in 
Canada where he started his career 
with theoretical and observational re¬ 
search. Later, he focused on instrumen¬ 
tation, in particular the construction 
of large telescopes. On the Canadian 
side, Odgers was instrumental in the 
construction of the CFHT on Mauna 
Kea. 
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astrophysical community: The circumstances under which multi¬ 
ple solutions occurred were either stellar-physically unacceptable 
or they were so restrictive that their realization in nature seemed 
unlikely A further handicap of the Odgers paper was that it was 
published and circulated as an observatory report only and as such 
must have had a diminished audience. Last but not least, the re¬ 
port's content is very formal and likely was too tough to digest for 
many in the astronomical community. 

In the second edition of his textbook, Aufbau und Enwicklung der 
Sterne, Vogt ( [1957^ came back to the Eindeutigkeitssatz der Theorie 
des Sternaufbaus and devoted a whole chapter to is. As in earlier 
publications, Vogt did not adopt a particularly mathematical point 
of view to discuss the problem. Instead, he merely insisted that 
the structure equations admit of unique solutions as long as the 
material properties (such equation of state, opacity, or nuclear 
energy generahon) are well defined functions of thermod3mamic 
state variables. As means to destroy the one-parameter families 
(characterized by different chemical composition) of evolutionary 
tracks on the HR plane, Vogt contemplated physical effects such as 
rotation, electromagnetic braking, or tidal effects in binaries. 

Despite the lack of evidence of any impact in stellar astrophysics 
of the report of Odgers ( [1957^ his paper marks the begin of the 
era in which multiple solutions began to pop up in the ever more 


detailed stellar-model computations. For example, Cox & Salpeter 


( [1964^ calculated low-mass pure helium star models in complete 
equilibrium. Below a critical mass, Mmin/ helium stars cannot 
maintain steady helium burning. Around this Mmin ~ O.SOSM©, 
a low- and a high-density solution with the same total mass and 
with identical composition were revealed. Not much later, for 
pure carbon stars too, double solutions for equal-mass model stars 
were encountered ( Deinzer & SalpeterjigB^ . For carbon stars, the 
minimum mass is larger than that of the helium stars. Neither in 
Cox & Salpeter ([igBqll nor in Deinzer & Salpeter is there any 


indicahon that the authors connected the double solutions with a 
failure of the VR theorem. Bodenheimer ( |i966[ l on the other hand 
questioned the validity of the VR theorem upon realizing that his 
pre - main-sequence model stars all converged essenhally to the 
same evolutionary locus along the Hayashi line, independent of 
the initial conditions he prescribed for his model sequences. Even 
though Bodenheimer was the most attentive author back then, his 
models do not serve as counterexamples to the VR theorem because 
they are not in complete equilibrium as they need to be for the 
orginal VR theorem to be applicable. 

A next higher level of complexity in stellar modeling was 
reached with composite models that consisted of a core and of 
a grafted envelope, both in complete equilibrium but both with a 
differing chemical composihon. Adopting the mass of the core as 
the control parameter allows to study the dependence of the physi¬ 
cal properties of a series of star models under conhnuous variation 


At the time of this writing, an ADS 
query returned 6 citations to the report; 
all between 1972 and 1978. 


No neutrino losses: Mmin ~ O.7M0, 
including neutrinos: Mmin ~ 0 . 82 Mo 
for steady carbon burning. 
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of the control parameter, such model sequences are known as linear 
series. To investigate the onset of the Schonberg-Chandrasekhar 
instability, [Gabriel & Ledoilx ( 1967^ chose the relative mass, 
of the inert helium core below a hydrogen-burning shell and a 
hydrogen-rich envelope as the control parameter of their linear 
series of model stars with constant total mass . Gabriel and Ledoux 
concluded that the instability develops at a turning point of their 
linear series. In the neighborhood of this turning point, the stellar 
structure equations were observed to admit of double solutions at 
constant (jHe- Investigations of the stability of these double solu¬ 
tions revealed then that one branch was secularly unstable. Even 
though double-solutions for the same stellar mass and the same 
chemical-composition profile were encountered, the result was 
not yet discussed in the context of the VR theorem. The situa¬ 
tion changed when Gabriel & Noels-Grotsch ( 1968^ studied pure 
carbon stars in the neighborhood of the respective Mmin/ again 
they found that only one branch, the low-degeneracy one, of the 
double-soluhon region was secularly stable. Eventually, the authors 
concluded that turning points of linear series signal violations of 
the classical VR theorem. Resorting to a more restricted formula¬ 
tion, [Gabriel & Noels-Grotsch] ( |i968^ tried to save the VR theorem 
by adding the aspect of secular stability: [...] For a given mass and 
chemical composition there exists only one secularly stable configuration. 
[...] 


Gonsul ting [Cox & Giuli , who published a comprehensive 

textbook which details knowledge and understanding of structure 
and evolution of simple single stars by the mid 1960s, one finds 
it to offer indeed an appropriate endpoint to the first part of this 
review of the history of the VR theorem. In PSS, a whole - albeit 
short - chapter is devoted to the VR theorem; that choice met criti¬ 
cism already early on by one of the reviewers of the books ([Sweet 


1969) and it likely sheds light on the importance attributed to the 
VR theorem at that hme. Be it as it may, Ghapter 18 of PSS is very 
useful here because it offers a glimpse at the perception of the VR 
theorem in the mid 1960s. Early on in the discussion, the authors 
emphasized that the VR theorem is not a theorem in a strict sense 
because cases of multiple solutions had been encountered and that 
a watertight proof had never been put forth. Nonetheless, Gox and 
Giuli could not resist the temptation to give a kind of a plausibility- 
'proof' of the VR theorem, following the line of thought already 
present in RusseU ( |i93ia| |. The system of equations that entered the 
proof remained those of a stellar configuration in complete equilib¬ 
rium so that the problem reduced to system of ordinary differential 
equation. The separated boundary conditions were introduced 
and it was argued, without going into any mathematical detail, 
that the implied algebraic constraints of the boundary conditions, 
being of lower dimensionality than the dimension of the solution 
space, ensure [...] under ordinary conditions [...] unique solutions. 
More originally, Gox and Giuli offered also a physical interpreta- 


N.B. Increasing the magnitude of (jHe 
can be understood as an emulation of 
stellar evolution through a sequence of 
equilibrium states. 


5 The monograph. Principles of Stel¬ 
lar Structure, is referred to as PSS 
subsequently. 

“ Simple in the sense of radial sym¬ 
metry of the star models, devoid of 
rotation and magnetic fields. 


The known double-solutions were 
either rather abstract and of little 
relevance to the Universe or one of 
the double solutions was secularly 
unstable and would not survive for 
long in nature. 
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tion of the VR theorem, they resorted to an order-of-magnitude 
discussion of the stellar structure problem (PSS, chapter 18.2) and 
thereby collapsed the differential equahons to a set of algebraic 
relations. They showed that if pressure, density, temperature, and 
radiated luminosity were prescribed, all emerging relations can be 
expressed as functions of mass, radius, and chemical composihon. 
Upon prescribing additionally also thermal equilibrium, the radius 
dependence of the set of algebraic equations can be eliminated so 
that eventually the order-of-magnitude approximations of the phys¬ 
ical quantities of a star in hydrostatic and thermal equilibrium are 
found to depend on mass and chemical composition only. Although 
intuitively attractive, the method applied only to a model star in a 
coarse integral sense and failed to be mathematically rigorous (cf. 
Appendix B). 

Astrophysics is not mathematics; in the latter, the VR theorem 
would have met its fate once one counterexample popped up - in¬ 
dependent of how academic it were. The astrophysical commimity, 
on the other hand, put up with the dilemma of the VR theorem and 
its counterexamples. After all, in astronomically relevant cases it 
seemed to remain predictive and explanations sounded plausible. 
Nonetheless, the beginning era of ever faster and easier accessible 
electronic computers at the end of the 1960s allowed to compute 
physically complex models in large numbers and the rapidly grow¬ 
ing repository of stellar models had intriguing challenges in store. 

The forthcoming second part of this marginalia on the history of 
the VR theorem will focus on the developments in the 1970s when 
a few astrophysicists set out to look more closely into matter of 
the VR theorem after ever more complicated star models could be 
computed and some of them exposed violahons of the VR theorem. 
The field benefited from fresh insights imported by people who 
applied to stellar astrophysics heavier mathematical machinery 
than usual. 


Appendix A: Eddington's mass - luminosity relation 

A correlation of stars' masses and luminosihes was hinted at as 
early as 1911; it was first mentioned in a clause of a paper on 
the motion of the stars in the Galaxy ( |Halm|i9ii[ |. Analyzing 
an appropriately chosen sample of 14 binary stars, [Hertzsprung 
( |i923^ was able to report quantitahvely on a relation of mass and 
brightness of his sample stars, finding clear evidence that more 
massive ones were consistently brighter than the less massive 
brethren. It was Eddington, however, who pushed the issue of the 
mass - luminosity relahon (M — L relahon in the following) further 
because he relied on it as an observational foundahon on which he 
could rest his mathemahcal modeling of the internal structure of 
the stars (|Eddington| 1924I1 . 
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Figure [^displays the same data as were shown in Fig.i of Ed- 


[dington] ( |i924) . The different markers in the plot idenhfy different 
classes of calibrated stars; filled circles show the pb stars, the pri¬ 
mary binary stars, triangles are the s b, the secondary binaries, 
squares show the eb, the eclipsing binaries, and finally, the asterisks 
stand for cep, the pulsahng variables. Eddington referred to them 
collectively as Cepheids; this despite the fact that the star RR Lyr 
was also in the sample. The physical differences between Cepheids 
and RR Lyrae variables were not known at the time. 

Erom the present point of view, Eddington's adding pulsating 
stars, the cep group, to the graph to make the case of an M — L 
relation is a dubious undertaking. It was essentially pure luck that 
made the outcome to look so seemingly convincing. At the time of 
Eddington's article, no direct determinations of masses of Cepheids 
were available. Eor the stars to find their place in the figure, Ed¬ 
dington resorted to his theorehcal M — L relation and applied it ad 
hoc to the Cepheids too. The computed stellar mass was then iter¬ 
ated imtil the pulsahon period of the modeled Cepheid eventually 
converged to the observed period. Therefore, the pulsating stars 
in Eig.|^are no fundamental indicators of the existence an M — L 
law but the result of an already plugged-in M — L relahon coupled 
with the pulsation theory of |Eddrngton| ( [i9i9| |. In contrast, the 
derivation of the masses of the members of the binary-star sample 
relies on Newton's laws acting in a 1 / r-gravitational potential only. 
Therefore, only the binary stars serve as legitimate indicators of the 
correlation between mass and luminosity. 

An inspection of the right panel of Eig.|^ makes clear that correla¬ 
tion of M — L relation of the binary stars persists also with modern 
physical calibrations of the binary stars, although the scattering 
increases at very low and very high luminosities. Eurthermore, 
the relation based on modern data has a steeper slope plus a few 
outliers in the graph; the Cepheids in particular scatter. Concen- 


Figure 3: Mass-luminosity relation 
as used in Eddington's studies in the 
early 1920s. The left panel contains 
the stars with the calibration known 
to |Eddington| 11924} . The right panel 
shows the relation when using current 
calibrations of the same objects. (Data 
harvested via SIMBAD and ADS). 
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trating on the arguably small number of Cepheids with the equally 
small mass spread, the modern data do not really call for the same 
M — L relation as it is suggested by the mostly main-sequence binary 
stars. Nonetheless, the overall relation is still impressively tight, 
particularly when accounting for the fact that very different kinds 
(evolutionary stages) of stars convene in the graph. 

Because the stars spend most of their lifetime buring hydrogen 
most stars observed in the sky are therefore likely in their main- 
sequence phase. Along the main sequence, the stars' luminosities 
grow with increasing mass. One representative form of an empir¬ 


ically calibrated main-sequence M — L relation is e.g. from Smith 


logL/Le =4.O-logM/M0 for M/M© > 0 . 43 . 


The binary-star data entering Eddington's M — L relation fit the 
above relation quite well so that we can assume that the respective 
stars are indeed likely in their main-sequence phase of their life. 

The Cepheids, on the other hand, are as we know today radi¬ 
ally pulsating intermediate-mass supergiants. Under favorable 
circumstances intermediate-mass stars loop across the HR diagram 
during their central helium-burning stage and some of these loop¬ 
ing stars migrate through the classical instability strip to become 
then observable as Cepheids. From stellar-evolution modeling we 
learned that close to the instability strip the luminosities of the 
same branches of these blue loops tend to be ordered in mass; i.e. 
Cepheids of different mass but a comparable evolutionary stage 
tend to obey an M — L relation too. The 2 "*^ crossing of the insta¬ 
bility strip during the early core helium-burning phase is usually 
the slowest and therefore the favored one to observe Cepheids in. 
Adopting hence this second crossing as the relevant one here, a 
fit to the computed intersections of evolutionary tracks with the 
instability strip ( [Chiosi et al.|i993| reads as 

logL/L0 = 3 . 57 -log M/M0-h 0 . 54 . 

Interestingly enough, the slopes of the Cepheids' and the main- 
sequence stars' M — L relation happen to be quite similar. From 
all we know, this is an accident of nature. The vertical displacement, 
A log(L/L0), of the two relations say at M* = 5 M© is only 0.23 so 
that the composite nature of the observed M — L relation in Fig. 
is hardly discernible, in particular in the presence of unavoidable 
scattering of observational data. 

Apart from some increased scattering in the modern version of 
Eddington's M — L relation, introduced by the pulsating variables, 
one obvious disagreement is apparent in comparison with the orig¬ 
inal one: The data point of the star RR Lyr lies far off the general 
trend. This is no surprise: Today we know that the family of RR 
Lyr variables is made up of low-mass (w O.6M0) population 11 stars 
living on the horizontal branch. Rather than following an M — L 
relation of the above kind, the spread in mass and in luminosity is 
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small SO that the class of RR Lyrae variables would form kind of a 
clump around the isolated prototype RR Lyr in Fig. (right panel). 


Appendix B: The nature of the equations 


To properly state what astrophysicists mean if they talk stars on 
the theory level, the set of the governing equations and assump¬ 
tions are laid out in the following. In each case of the following 
collection of formulae, the first line contains, as the starting point, 
the general fluid-d 3 mamical equations in their Lagrangian form. 
We follow mostly the notation of Mihalas & Mihalas-Weibel ( 1984 ^ 
(only if not self-evident, deviations therefrom are explained). The 
second line specializes then on spherical symmetry^ ^ and on mass 
as the independent variable - choices usually adopted in stellar 
structure/evolution modeling : 


Dtp = -p (V ■ p) , 

( 1 ) 

°'”"“47rrV 

( 2 ) 


AO = AnpG, 

( 3 ) 

Gm 

(4) 

S= 


In spherical symmetry, eq. |^can be integrated once; defining g = 
—d,xl> leads then to eq.|^ 


pDtV = /-h V ■ T, 

(5) 

Gm 1 0 

DmP - . 4 . 2 ■ 

4:71 4:71 

(6) 


The stress tensor T is made up of the components of viscous stress 
and it harbors on its diagonal the isotropic hydrostatic pressure. 
Neglecting viscosity turns Cauchy's equation into Euler's equation. 
Physically, the stress tensor is relatively imimportant in stellar 
matter; for numerical reasons (in 'hydrodynamicaT computations) 
an artificial stress tensor might become instrumental to smear out 
shocks and to stabilize the computational scheme. 

The change of the specific heat content of the stellar matter, q, 
remodeled with the help of the 1®* law of thermodynamics, can be 
written as: 

~ ^visc V ■ F -)- p Snuc / (7) 

Dfp = —D^F “he. (8) 

As in the momentum equation, the dissipation term - here the vis¬ 
cous energy dissipation function Yyisc “ is physically not important 
in stellar interiors except if sharp fronts develop and then in con¬ 
nection with the necessity to stabilize the numerical treatment. In 


pDtq = p 


Dte + PDt [ - 
P 


” Spherical symmetry is appropriate 
for non-rotating, non-magnetic star. 
Even though it seems intuitively 
obvious that a static self-gravitating 
fluid configuration assumes the form 
of a sphere, the proof that the sphere 
is the only equilibrium figure came 
relatively late, see e.g. [Poincare & 
|Dreyfus|(i902[ with a proof which 
relied on Lyapunov's master thesis of 
1884 or consult [Carleman|(i9i9) who 
used a more geometrical ansatz. 

Continuity equation 
Mass equation 


Poisson equation (elliptic equation) 


Cauchy's equation (hyperbolic equa¬ 
tion) 


Gas energy equation 
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canonical quasi-static stellar-evolution computations, Yyisc can be 
neglected. 

The energy flux is denoted by F. Astronomers prefer to work 
with the local radially streaming luminosity, which is defined as 
L = inr^Fr, with Fr being the radial component of the local energy 
flux. Furthermore, the energy source, s, is attributed to nuclear 
burning in the star, therefore it is subscripted with 'nuc'. In the 
stellar astrophysical form of the equation, as shown on the second 
line, the nuclear energy input rate, measured as energy per unit 
time and unit mass, is referred to as e. It is important to keep in 
mind that energy gain by nuclear burning contributes positively 
and possible energy loss by neutrino production is to be subtracted 
because neutrinos do not contribute to the heat content of the 
stellar matter. 

To get a handle on the physics of the energy flux, stellar astro¬ 
physicists adopt Fourier's law to model the flow of photons as 
a diffusion process driven by the spatial temperature gradient. 
Even energy transport by material motion can be appropriately 
accommodated. 


Extreme conditions such as encoun¬ 
tered e.g. during a stellar core collapse 
can trap even neutrinos and hence 
require then a careful treatment of the 
energy budget including the neutrinos. 


F=-K-VT, (9) 

F>mT = - i^Tir^p ■ Vo- (lo) 

The quantity K denotes the coefficient of thermal conductivity; 
in radiative regions it can be written as K = acT^/ (3Kp) with 
k(p, T,x) being the Rosseland opacity. The quantity Vq measures 
the temperature stratification: 


Fourier's law for the flux (parabolic 
equation) 


din T 


Vrad radiative region, 
Vc convective region. 


with Vrad = 3kLP / {16nacGmT^) being a purely local function of 
stellar quantities. In the easiest case of a stellar-convection descrip¬ 
tion, such as in elementary mixing-length models for example, also 
Vc is a function of local variables alone. More elaborate treatments 
of convection can, however, introduce non-local contributions. 

Nuclear burning is the source of a star's evolution; the resulting 
spatio-temporal change of nuclear species X, in a star can formally 
be written as 

DfX,- = Qi -Si + V (udVX,) . (11) 

Apart from the source-, Q/, and the sink-term, S,, both determined 
by the type and complexity of the nuclear burning network, nuclear 
species can be smeared out spatially by a multitude of physical pro¬ 
cesses (such as convection, thermohaline mixing, semi-convection, 
settling, levitation); these transport processes are hidden away in 
a diffusion-type term in the equation with the particular physi¬ 
cal process manifesting itself in the specification of the diffusion 
coefficient up). In the majority of the numerical realizations, the 


Xj be the mass fraction of nuclear 
species i; E; X/ = 1. 

i e [l, ..., Nspec], with Nspec the 

number of species accounted for in the 
stellar matter. 
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computation of the nuclear evolution is decoupled from the stellar 
structure problemd^ We presume that at each epoch t, the vec¬ 
tor X = (x’l {m, t),..., {m, t)'^ is known via some suitable 

computational procedure. 

Boundary conditions for the stellar structure equahons are dis¬ 
tributed ones with the natural choices in the center: r = 0 and L = 0 
at m = 0 . The surface is, by its very stellar nature, ill-defined and 
requires suitably chosen physical approximations: Traditionally 
popular is the assumphon of thermal equilibrium of radiation 
and matter fields at the photosphere leading to: L = inr^aT^ at 
m = M*. The radius at the photosphere is then set equal to star's 
radius r = R^, and the temperature at the photosphere corresponds 
to the so called effective temperature T'eff- The second outer boimd- 
ary condition, a mechanical one, determines e.g. the pressure at the 
photosphere: 

P=f{p, T,k{p, T,x)), 

with some suitable function /, which approximates the type of 
atmosphere which exerts its pressure, P, on the photosphere. For 
simplicity's sake, and likely sufficient for pure mathemafical consid- 
erafions, it is sufficient to assume some ad hoc constant pressure at 
the photosphere: 


" For a modern approach coupling 
nuclear with stellar-structure evolution, 
see Appendix B of|Faxton et al.|(20i3) 


P — f^phot — const. <SC Tcenter • 


Finally, initial data that specify the state of the star's structure are 
required to get a model sequence started in time. Typically, such 
a hme evolution is initialized with some simplified star model in 
hydrostatic and if possible also in thermal equilibrium. Both as¬ 
sumptions ensine that pressure and temperature are continuously 
differentiable in space. Density, on the other hand, can have dis¬ 
continuities, depending on the spatial structure of the composition 
vector; luminosity will also react accordingly. Think, for example, 
of an initial model with a pure helium core and a pure hydrogen 
envelope: Across the H/He interface pressure and tempera hire are 
continuous whereas density and, at sufficiently high temperature, 
also luminosity develop finite jumps. 


Where astronomers were apparently too light-hearted in 'proof¬ 
ing' the VR theorem, mathematicians, on the other hand, in par¬ 
ticular those rising a warning finger, such as Kurfh ( [i953[ l, were 
essentially absent. The asterophobia of the mathematicians is likely 
caused by the fact that the equations that model the structure of 
the stars cannot be pigeon-holed: Depending on the specific sim¬ 
plifications introduced to the system of structure and evolution 
equations, they can change the mathematical character so that dif¬ 
ferent mathematical tools must be applied to the formal study the 
problem; on the other hand also different numerical methods must 
be implemented to tackle the computational problem. The simplest 
approach to model stars, which likely has canalized early think¬ 
ing of proofing the VR theorem, namely the one used to compute 
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polytropes is shortly touched upon in the following. 


Separating mechanical and thermal parts of the stellar 
structure problem was the first ansatz to come to grips with under¬ 
standing the interior conditions of stars. To establish the necessary 
barotropic conditions, frequently a polytropic^^ relation was postu¬ 
lated. In the static case, i.e. in absence of a velocity field, eqs. 2, 4, 
and 6 morph into the venerable Lane-Ritter-Emden equation. In the 
formative years of fheorefical asfrophysics, fhis equation was solved 
as an initial-value problem (IVP): The computation started in the 
regularly singular center of fhe model wifh prescribed values of fhe 
dependenf variable and ifs derivative. The infegrafion was followed 
out to the first root of the dependent variable; its location was then 
adopted as a measure of the radius of the model star. Looking at 
the problem as an IVP ensured existence and uniqueness of fhe 
solution by the sufficiently smooth character of the right-hand side 
of fhe ODE via Picard-Lindelof's fheorem. 

More generally, in particular with a non-vanishing velocity 
term, the original mechanical fluid-dynamical equations (eqs. 

1, 3, and 5, closed with a polytropic relation between p and P) 
constitute the so-called Euler-Poisson problem. Assuming a compact 
support for the density and hence for the whole problem, i.e. p > 0 
obtains for a finite volume only and defining the outer boundary 
of the gravitationally-bound fluid sphere by p(R) = 0 bring about 
considerable mathematical complications; a substantial body of 
literature exists on uniqueness and evolution of such boundary- 
value problems (BVPs). Eor some recent advances, consult e.g. 
Deng & Guo ([2003[| who proofed uniqueness theorems for the 


static case, i.e. the BVP which must be solved for the structure of 
polyfropic spheres. 

Even though the proof of the VR theorem in Russell ( [i93ia[ | 

(and all later repetitions thereof) refers to the stellar-structure 
problem as a BVP, the presentation of how fhe equafions are solved 
is reminiscent of the direct integration of the Lane-Ritter-Emden 
equation. In other words, from reading Russell ( [i93ia[ | one comes 
away wifh the impression that existence and uniqueness properties 
were implicitly influenced by the IVP experience gained with 
polytropes. However, existence and uniqueness statements for 
higher-order BVPs are mathematically formidable; we superficially 
referred fo the easiest case, that of the Euler-Poisson system, just 
before. 


’3 The relation P a with n 

being the polytropic index, constitutes 
a stratification relation rather than a 
state relation within a prescribed fluid 
element. 


Properties of the equations of more realistic approximations to 
the stars' structure and their evolution will be a topic in the second 
part of fhis essay on fhe history of fhe Vogf-Russell fheorem. 
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